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Abstract 

SOMIM is an open-source code that implements a Search for 
Optimal Measurements by using an Iterative Method. For a 
given set of statistical operators, SOMIM finds the POVM that 
maximizes the accessed information, and thus determines the 
accessible information and one or all of the POVMs that retrieve 
it. The maximization procedure is a steepest-ascent method that 
follows the gradient in the POVM space, and uses conjugate 
gradients for speed-up. 
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1 License Agreement 

SOMIM is an open-source code that, for a given set of statistical opera- 
tors, implements a Search for Optimal Measurements by using an Iterative 
Method. © 2007 K.L. Lee, W.K. Chua, S.Y. Looi and B.-G. Englert. 

SOMIM is free software: You can redistribute it and/or modify it under 
the terms of the GNU General Public License Version 3 as published by the 
Free Software Foundation. 

SOMIM is distributed in the hope that it wih be useful, but WITH- 
OUT ANY WARRANTY; without even the implied warranty of MER- 
CHANTABILITY or FITNESS FOR PARTICULAR PURPOSE. See the 
GNU General Public License at http : / /www . gnu . org/licenses/, for de- 
tails. 

2 What can SOMIM be used for? 

Given a set of statistical operators pj, SOMIM calculates the accessible 
information (AI) associated with these statistical operators, and finds the 
POVM that retrieves the AI, or rather: it finds one of the optimal POVMs, 
as the optimum need not be unique. Repeated runs of SOMIM for the 
same input data but different random seeds may yield alternative optimal 
POVMs. 

The calculation is performed using the steepest-ascent iteration intro- 
duced in Ref. [Ij; see also Section 11.5 in Ref. [2]. The implementation in 
SOMIM makes use of conjugate gradients and the Golden-Section search 
method. 
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3 Download and compile (optional) 

The complete set of files, including this manual, are available at the SOMIM 



site: |http : / /theory . quantumlah . org/pro j ect/SOMIM/ Download http 



//theory. quantumlah.org/project/SOMIM/all. zip, if you want to have 
the complete collection of files. Just this manual is fetched from http : ' 
//theory . quantumlah . org/pro j ect/SOMIM/SDMIMmanu al . pdf The Win- 



dows executable file for SOMIM can be downloaded from http: //theory 



quantumlah.org/project/SOMIM/somim.zip If you intend to modify the 



code, you could download the source files from jhttp : //theory . quantumlalT 
[org/pro j ect/SOMIM/ source . zip, 

The program is written in C-|--|- and the graphic user interface (GUI) 
is implemented using wx Widgets (http://www.wxwidgets.org/). Here are 
the instructions for compiling SOMIM: 



1. Install wx Widgets from http : //www . wxwidgets . org/ downloads/] 



2. If you are working in Windows, you need to install MinGW (http:// 



www.mingw.org/download.shtml) and MYSY (http: //www. mingw, 
org/msys . shtml) as well. 



3. When wxWidgets and MinGW are configured, you can compile SOMIM 
by executing gH — |- MI.cpp 'wx-config libs' 'wx-config cxxflags' 
-o YourProgramName in MSYS shell. 

4. If you face problems running the program in a Linux environment, try 
export LDJLIBRARY PATH=/usr/local/lib 

5. The executable file is compiled under Windows XP Service Pack 2, 
with wxWidgets 2.8.6, MinGW 5.1.3, MSYS 1.0. 



4 How to use the program? 

The GUI of SOMIM is shown in Fig.[Tj In the first box labeled as "Parameter 
Settings", J is the number of statistical operators to be input. The current 
maximum possible value is J = 30. You can change the maximum value by 
modifying the source code. Parameter K is the initial number of elements of 
POVM, with the largest possible value being K = 30. The third field is the 
dimension N of the statistical operators which has = 30 as the highest 
possible value; that is: the J different statistical operators are represented by 
N X N matrices. The fourth field is the chance of using the steepest- ascent 
gradient method to perform maximization in an iteration; this parameter 
controls the relative frequency of using the direct or the conjugate gradient. 
The fifth field gives the tolerance in the accessible information, the stopping 
criteria for the computation; The calculation stops when the difference in ac- 
cessible information between the current iteration and the previous iteration 
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4 HOW TO USE THE PROGRAM? 



File Calc About 
test3.txt 

Parameter Settings- 



Number of possible quantum stateSj J = 
Number of POVM, K = 
Dimension of statistical operator rho J - 
Percentage chance to calculate with direct gradient = 
Tolerance in Mutual Information = 

Output file = output.dat 



X 4 



1° 

0.000001 



Input the Matrix Representation of Quantum States - 
Input the matrix representation of rho 4 



-Display Real or Imaginary 

C* Real C Imaginary 



1^ 0.1938 


0,0348 


0.0348 


0.0348 


0.0348 


0,0187 


0.0000 


0.0000 


0.0348 


0,0000 


0.0187 


0.0000 


0.0348 


0,0000 


0.0000 


0.0187. 



pDisplay Optimal POVMs - 



Display the optimal POVM PI 1 1 



■Display Real or Imaginary 

!'• Real Imaginary 



1 0.2501 


-0,2500 0.2501 0.2500 


-0.2500 


0,2499 


-0.2500 


-0.2499 


0.2501 


-0,2500 


0.2501 


0.2500 


0.2500 


-0,2499 


0.2500 


0.2498 



-Display the maximal Mutual Information - 



Maximal Mutual Information^ 0.232544 



Calculate MI 



Reset 



Figure 1: The graphic user interface (GUI) of the program. 
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is less than half of the sum multiplied by the tolerance plus 10~ , i.e. when 
2.0 X (current — previous) < tolerance x (current + previous) + 1.0 x 10~^^. 
The sixth field is the file name of the output file. The output file will be 
located in the same directory as the program. 

The next two boxes display the statistical operators {pj}j=i,...,j and the 
calculated elements of the optimal POVM {Ilk}k=i,...,K- The spin buttons 
are used to switch between the various pjs and II^s, respectively, while the 
small box beside the spin button is used to choose to display the real or 
imaginary part of the chosen pj or U^. 

The maximum accessible information for the given set {pj} will be dis- 
played in the last box after the "Calculate MI" button is pressed. All values 
will be reset to default when the "Reset" button is pressed. 

Important note: The matrices for the pjS must have the correct dimen- 
sion; they have to be hermitian with nonnegative eigenvalues; their traces 

J 

are their statistical weights, which must add to unity: ^tr{pj} = 1 . 

i=i 

5 How to import data? 

Data can be imported into SOMIM using a text file that is possibly generated 
by another program. An example is shown in Fig. [2} When importing, the 
numbers after the equal signs will be read into SOMIM. The first line gives 
the dimension N of the statistical operators. The second line gives the 
number J of statistical operators while the third line gives the number K 
of elements of the POVMs that SOMIM should start calculating with. 



r — 

A.test3.txt - Notepad 






_ n 


1 

X 


File Edit Format View 


Help 








N= 4 
3= 4 
K= 7 

rhol = {{0.19375, 
rho2 = {{0.19375, 
rho3 = {{0.19375, 
rho4 = {{0.19375, 


0.0347985, - 
-0. 0347985, 
-0. 0347985, 
0.0347985, 


0. 0347985, 
.0347985, - 
-0. 0347985, 
.0347985, 


-0.0347985}, {0.034 
0.0347985}, {-0.034 
0.0347985}, {-0.03. 
.0347985}, {0.03479 




.<J.. . . .1111 






3 





Figure 2: Example of an import file. 

Subsequent lines give the input matrices for the statistical operators. 
Each line will give only one operator. For an operator represented in matrix 
form as 

0.1 0.3 + 0.5i \ , . 

0.3-0.5i 0.6 j ' ^ ' 
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7 CONTACT INFORMATION 



the input data should be formatted as {{0.1, 0.3 + 0.51}, {0.3 - 0.51, 0.6}}. 

Complex numbers are entered as RealPart + ImaginaryPart I, as illus- 
trated by —3.1 — 4.51. Please note that the complex unit i must be entered 
in upper case (I) and it must be at the end of the entry. 

6 Meaning of output data 

A typical output file looks like Fig. |3j The first six lines give the follow- 
ing information: the number J of statistical operators, the initial number 
K of elements in POVM, the dimension of the statistical operators, the 
probability of maximizing information using steepest gradient method com- 
pared to conjugate gradient method, the tolerance in the calculated mutual 
information, and the seed for the random number generator. 

The next block of lines gives the mutual information at the end of each 
iteration, 

J K 

new_MI = V V Pjk log2 

with 

K J 

Pjk = tiiPjUk} , Pj. = Pjk , p-k = Yl ■ 

k=l 3=1 

In the example shown in Fig. |3] 27 iterations have been performed with the 
final accessible information being AI = 0.232 544. 

The subsequent two blocks of lines give the J statistical operators and 
the K elements of the optimal POVM that correspond to the accessible 
information calculated in the final round of iteration. Each statistical oper- 
ator /POVM element is given in a single line in matrix form, as explained in 
Section IH 

Among the, K POVM elements, if any two elements, let's say 11^^ and 
11^2, give equivalent probabilities, i.e. PjkiP-k2 = P-kiPjk2 for J; then 
these two POVM elements are replaced with one new POVM element, n^i -|- 
11^2, such that the new optimal POVM contains only K — 1 elements. The 
last block of data in the output file gives the POVM after this elimination 
process, i.e. the POVM is the optimal POVM with the least number of 
elements. 

7 Contact information 

Please send your comments, suggestions, or bug reports to the following 
email account: somim@quantumlah.org 
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W output, ddt - Notepad 



File Edit Format View Help 



rho_j = 4 
POWl = 7 
of rho_j = 4 
direct gradient = 0/100 
in Mutual Information = le-006 
seed for number generator = 1197369744 
1, new_MI=0. 00926005 
new_MI=0. 119213 
new_MI=0. 172292 
new_HI=0. 18576 
new_MI=0. 198619 
new_MI=0. 220146 
new_HI=0. 229224 
new_MI=0. 23038 
new_MI=0. 230921 



10, 


new_MI 


=0 


231371 


11, 


new_MI 


=0 


232025 


12, 


new_Ml 


=0 


232098 


13, 


new_Ml 


=0 


232189 


14, 


new_Hl 


=0 


232355 


15, 


new_J^1I 


=0 


232448 


16, 


new_Ml 


=0 


232471 


17, 


new_Ml 


=0 


232497 


18, 


new_MI 


=0 


232512 


19, 


new_Ml 


=0 


232518 


20, 


new_Ml 


=0 


232526 


21, 


new_Hl 


=0 


232537 


22, 


new_j^ii 


=0 


23254 


23, 


newjil 


=0 


232541 


24, 


new_Ml 


=0 


232542 


25, 


new_MI 


=0 


232543 


26, 


new_Ml 


=0 


232543 


27, 


new_Ml 


=0 


232544 



Number of 

Number of 

Dimensi on 

chance of 

Tol erance 

starti ng 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

counter 



rhol={{0.1937 5, 0. 034798 5,-0. 034798 5 , -0 . 0347985 }, {0. 0347985 , 0. 01875 , 001 , O+OI }, {-0. 0347 
rho2={{0.1937 5, -0. 034798 5, 0. 034798 5 , -0 . 0347985 }, {-0. 034798 5, 0. 01875 , OOl , OOl }, {0. 034 79 
rho3={{0.1937 5, -0. 034798 5, -0. 0347985, 0. 0347985}, {-0. 034798 5, 0. 01875 , O+OI , O+OI }, -[-0. 03. 
rho4={{0.1937 5, 0. 034798 5, 0. 034798 5, 0. 034798 5}, {0. 034 7985, . 0137 5 , O+OI , 001 }, {0. 0347985 
POVW1={{0. 2 50265, -0. 2 50027+3. 1958e-005 I, -0. 2 50115+0. 00015 96531 , -0. 2 5 0121+0. 000133 5 941 

POWI2={{0. 241149, 0. 241192-2. 98902 e-005 1 , -0. 241318+0. 0001688811 , 0. 241158-0 . 00024074 81 } 

POVm3={{0. 209 594, -0. 209 569-3 . 97619e-00 51 , 0. 2094 98-0. 000185 2621 , 0. 20949+6. 0026e-005I }, 

POVW4={{0. 249694, 0. 249899+4. 24432 6-0051, 0. 249914-0. 0001776381,-0. 24 9878+7. 52982e-00 5I 

POWI5={{0. 0405497, -0. 04 04033+ 5 . 10944e-0051 , 0. 04 5864 -2 . 70676e-0051 , 0. 040597-0. 0001181. 

POWI6={{0. 000336082, 0. 000346802+1 . 4 590 5e-00 51 , -0 . 0003 32003+2 . 0244 5e-006l , 0. 0003232 77+ 

POWI7={{0. 0084117, 0. 008 56016-7. 04 34e-00 51 , -0. 00823416+5 . 94103e-0051 , 0. 00843119+7. 7861 

After el imi ni nati ng the extra povms 

POWIl={{0. 2 50265, -0. 2 50027+3. 1958e-005 I, -0. 2 50115+0. 00015 96531 , -0. 2 5 0121+0. 000133 5 941 
POWI2={{0. 249897, 0. 2 50099-8. 57337e-005I , -0. 249884+0. 00023 03151 , 0. 24 9913-0. 000150771 }, 

■M III I >-i' 



^1 



Figure 3: Example of output file. 
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